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Detrended fluctuation analysis (DFA) and detrended moving average (DMA) are two scaling 
analysis methods designed to quantify correlations in noisy non-stationary signals. We systemati- 
cally study the performance of different variants of the DMA method when applied to artificially 
generated long-range power-law correlated signals with an a-priori known scaling exponent oq and 
compare them with the DFA method. We find that the scaling results obtained from different 
variants of the DMA method strongly depend on the type of the moving average filter. Further, 
we investigate the optimal scaling regime where the DFA and DMA methods accurately quantify 
the scaling exponent qq, and how this regime depends on the correlations in the signal. Finally, 
we develop a three-dimensional representation to determine how the stability of the scaling curves 
obtained from the DFA and DMA methods depends on the scale of analysis, the order of detrending, 
and the order of the moving average we use, as well as on the type of correlations in the signal. 
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component feedback interactions play a central role, ex- 
hibit complex self-similar fluctuations over a broad range 
of space and/or time scales. These fluctuating signals can 
be characterized by long-range power-law correlations. 
Due to nonlinear mechanisms controlling the underlying 
interactions, the output signals of complex systems are 
also typically non-stationary, characterized by embedded 
trends and heterogeneous segments (patches with differ- 
ent local statistical properties) J4, ^5, . Traditional 
methods such as power-spectrum and auto-correlation 
analysis Isof are not suitable for non-stationary 

signals. 

Recently, new methods have been developed to address 
the problem of accurately quantifying long-range corre- 
lations in non-stationary fluctuating signals: (a) the de- 
trended fluctuation analysis (DFA) 16, 23, 5 1| , and (b) 
the detrended moving average method (DMA) [52. [53. [54 
|5l[5i|. An advantage of the DFA method [I I 14511461 147; 
is that it can reliably quantify scaling features in the fluc- 
tuations by filtering out polynomial trends. However, 
trends may not necessarily be polynomial, and the DMA 
method was introduced to estimate correlation proper- 
ties of non-stationary signals without any assumptions on 
the type of trends, the probability distribution, or other 
characteristics of the underlying stochastic process. 

Here, we systematically compare the performance of 
the DFA and different variants of the DMA method. To 
this end we generate long-range power-law correlated sig- 



nals with an a-priori known correlation exponent ao us- 
ing the Fourier filtering method 57]. Tuning the value 
of the correlation exponent ao, we compare the scaling 
behavior obtained from the DFA and different variants 
of the DMA methods to determine: (1) how accurately 
these methods reproduce ao; (2) what are the limita- 
tions of the methods when applied to signals with small 
or large values of ao. Based on single realization as well 
as on ensemble averages of a large number of artificially 
generated signals, we also compare the best fitting range 
(i.e. the minimum and the maximum scales) over which 
the correlation exponent ao can be reliably estimated by 
the DFA and DMA methods. 

The outfine of this paper is as follows. In Sec. II, we 
review the DFA method and we introduce variants of 
the DMA method based on different types of moving av- 
erage filters. In Sec. HI we compare the performance 
of DFA and DMA on correlated and anti-correlated sig- 
nals. We also test and compare the stability of the scaling 
curves obtained by these methods by estimating the local 
scaling behavior within a given window of scales and for 
different scaling regions. In Sec. IV we summarize our 
results and discuss the advantages and disadvantages of 
the two methods. In Appendix I we consider higher or- 
der weighted detrended moving average methods, and in 
Appendix II we discuss moving average techniques in the 
frequency domain. 



II. 



METHODS 



A. Detrended Fluctuation Analysis 

The DFA method is a modified root-mean-square (rms) 
analysis of a random walk. Starting with a signal u{i), 
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where i = 1,...,N, and N is the length of the signal, 
the first step of the DFA method is to integrate u{i) and 
obtain 



where 
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is the mean. 

The integrated profile y{i) is then divided into boxes 
of equal length n. In each box n, we fit y{i) using a poly- 
nomial function ?/n(i), which represents the local trend 
in that box. When a different order of a polynomial fit 
is used, we have a different order DFA-£ (e.g., DFA-1 if 
e=l, DFA-2 if £ = 2, etc). 

Next, the integrated profile y(i) is detrended by sub- 
tracting the local trend ?/«(«) in each box of length n: 
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Finally, for each box n, the rms fluctuation for the inte- 
grated and detrended signal is calculated: 



F(>i) 
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The above calculation is then repeated for varied box 
length n to obtain the behavior of F{n) over a broad 
range of scales. For scale-invariant signals with power- 
law correlations, there is a power-law relationship be- 
tween the rms fluctuation function F{n) and the scale n: 



F{n) 
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Because power-laws are scale invariant, F(n) is also 
called the scaling function and the parameter a is the 
scaling exponent. The value of a represents the degree 
of the correlation in the signal: if a — 0.5, the signal 
is uncorrelated (white noise); if a > 0.5, the signal is 
correlated; if a < 0.5, the signal is anti-correlated. 



B. Detrended Moving Average Methods 

The DMA method is a new approach to quantify corre- 
lation properties in non-stationary signals with underly- 
ing trends [s^, El] . Moving average methods are widely 
used in fields such as chemical kinetics, biological pro- 
cesses, and finance [EE EE [59, 60, 61] to quantify signals 
where large high-frequency fluctuations may mask char- 
acteristic low-frequency patterns. Comparing each data 



point to the moving average, the DMA method deter- 
mines whether data follow the trend, and how deviations 
from the trend are correlated. 

Step 1: The first step of the DMA method is to detect 
trends in data employing a moving average. There are 
two important categories of moving average: (I) simple 
moving average and (II) weighted moving average. 

(I) Simple moving average. The simple moving 
average assigns equal weight to each data point in a win- 
dow of size n. The position to which the average of all 
weighted data points is assigned determines the relative 
contribution of the "past" and "future" points. In the fol- 
lowing we consider the backward and the centered mov- 
ing average. 

(a) Backward moving average. For a window of size n 
the simple backward moving average is defined as 



^ n—1 
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where y{i) is the integrated signal defined in Eq.(^. 
Here, the average of the signal data points within the 
window refers to the last datapoint covered by the win- 
dow. Thus, the operator y„ in Eq.© is "causal", i.e., 
the averaged value at each data point i depends only on 
the past n — 1 values of the signal. The backward moving 
average is however affected by a rather slow reaction to 
changes in the signal, due to a delay of length n/2 (half 
thewindow size) compared to the signal. 

(b) Centered moving average This is an alternative 
moving average method, where the average of the signal 
data points within a window of size n is placed at the 
center of the window. The moving average function is 
defined as 



yn{i) 
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where y(i) is the integrated signal defined in Eq.(Q and 
\x\ is the integer part of x. The function y„ defined in 
Eq.(|Zj) is not "causal", since the centered moving aver- 
age performs dynamic averaging of the signal by mixing 
data points lying to the left and to the right side of i. 
In practice, while the dynamical system under investiga- 
tion evolves with time i according to y(i), the output of 
Eq.Q mix past and future values of yii). However, this 
averaging procedure is more sensitive to changes in the 
signal without introducing delay in the moving average 
compared to the signal. 

(II) Weighted moving average. In dynamical sys- 
tems the most recent data points tend to refiect better 
the current state of the underlying "forces" . Thus, a filter 
that places more emphasis on the recent data values may 
be more useful in determining reversals of trends in data. 
A widely used filter is the exponentially weighted moving 
average, which we employ in our study. In the follow- 
ing we consider the backward and the centered weighted 
moving average. 
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(a) Backward moving average. The weighted backward 
moving average is defined as 



yn{i) = (1 - A)y(i) + Ay„(?; - 1), 



(8) 



where the parameter A = n/(n + 1), n is the window 
size, i — 2,3, ...,7V and yn(l) = 2/(1)- Expanding the 
term ?/„(z — 1) in Eq.(|Hl), we obtain a recursive relation of 
step one with previous data points weighted by increasing 
powers of A. Since A < 1, the contribution of the previous 
data points becomes exponentially small. The weighted 
backward moving average of higher order £ > 1 (WDMA- 
t) where ^ is the step size in the recursive Eq. ^ is defined 
in Appendix I. 

(b) Centered moving average. The weighted centered 
moving average is defined as 



(9) 



where is defined by Eq.®, and = (1— A)y(i) + 

\y^{i + 1), where i = iV - 1, - 2, 1 and ^^(A^) = 
y{N). The term is the weighted contribution of all 

data points to the right of i (from i + 1 to the end of the 
signal A^), and is the weighted contribution of all 

data points to the left of i. 

The exponentially weighted moving average reduces 
the correlation between the current data point at which 
the moving average window is positioned and the previ- 
ous and future points. 

Step 2: Once the moving average yn{i) is obtained, 
we next detrend the signal by subtracting the trend y„ 
from the integrated profile y[i) 



Cn{i) = y{i) - j/n(j)- 
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For the backward moving average, we then calculate the 
fluctuation for a window of size n as 
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For the centered moving average the fluctuation for a 
window of size n is calculated as 
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step 3: Repeating the calculation for different n, we 
obtain the fluctuation function F{n). A power law rela- 
tion between the fluctuation function F{n) and the scale 
n (see Eq.©) indicates a self-similar behavior. 

When the moving average y„ is calculated as in Eq.® , 
Eq.©, Eq.® and Eq.Q, we have the detrended moving 
method (DMA), the centered detrened moving average 
(CDMA), the weighted detrended moving average with 
order i (WDMA-^) and the weighted centered detrended 
moving average (WCDMA) respectively. 



(a) Anti-correlated signal: a„=0.2 
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(c) Correlated signal: a„=0.8 
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FIG. 1: A comparison of the scaling behavior obtained from 
the DMA, WDMA-1, DFA-0, and DFA-1 methods for arti- 
ficially generated power-law correlated signals with a scaling 
exponent ag. The length of the signals is A = 2^". Scaling 
curves F{n) versus, scale n for (a) an anti-correlated signal 
with Qo = 0.2, (b) an uncorrelated signal with ao = 0.5 and 
(c) a positively correlated signal with ao = 0.8. At small 
scales, all methods exhibit a weak crossover, which is more 
pronounced for anti-correlated signals. At large scales, the 
F{n) curves obtained from DMA, WDMA-1, and DFA-0 ex- 
hibit a clear crossover to a flat region for all signals, indepen- 
dent of the type of correlations. No such crossover is observed 
for the scaling curves obtained from the DFA-1 method, sug- 
gesting a more accurate estimate of the scaling exponent ao 
at large scales. 



III. ANALYSIS AND COMPARISON 

Using the modified Fourier filtering method [s^l, we 
generate uncorrelated, positively correlated, and anti- 
correlated signals u{i), where i = 1,2, N and A^ = 2^°, 
with a zero mean and unit standard deviation. By intro- 
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FIG. 2: A comparison of the performance of the different scal- 
ing methods (DMA, WDMA-1, DFA-0, DFA-1, and DFA-2) 
when applied to artificially generated signals with long-range 
power-law correlations. Here ao is the correlation exponent 
of the generated signals and a is the exponent value esti- 
mated using different methods. For all methods we obtain a 
by fitting the corresponding scaling curves F{n) in the range 
ne [f0^fO'']. Flat regions indicate the limitations of the meth- 
ods in accurately estimating the degree of correlations in the 
generated signals, eis the "output" exponent a remains un- 
changed when the "input" exponent ao is varied. 



ducinga designed power-law behavior in the Fourier spec- 
trum |45l Is?! , the method can efficiently generate signals 
with long-range power-law correlations characterized by 
an a-priori known correlation exponent ao ■ 



A. Detrended moving average method and DFA 

In this section we investigate the performance of the 
DMA and WDMA-1 methods when applied to signals 
with different type and degree of correlations, and com- 
pare them to the DFA method. Specifically, we compare 
the features of the scaling function F{n) obtained from 
the DMA and WDMA-1 methods with the DFA method, 
and how accurately these methods estimate the correla- 
tion properties of the artificially generated signals u(i). 
Ideally, the output scaling function F{n) should exhibit 
a power-law behavior over all scales n, characterized by 
a scaling exponent a which is identical to the given cor- 
relation exp onent oq of the artificial signals. Previous 
studies |4J, ]45, 46] show that the scaling behavior ob- 
tained from the DFA method depends on the scale n and 
the order £ of the polynomial fit when detrending the 
signal. We investigate if the results of the DMA and 
WDMA-1 method also have a similar dependence on the 
scale n. We also show how the scaling results depend on 
the order i when using WDMA-^ with i = 2,3,4,5 are 
applied to the signals (See Appendix I). 

To compare the performance of different methods, we 



first study the behavior of the scaling function F{n) ob- 
tained from DFA-0, DFA-1, DMA, and WDMA-1. In 
Fig. n we show the rms fiuctuation function F{n) ob- 
tained from the different methods for an anti-correlated 
signal with correlation exponent ao = 0.2, an uncorre- 
lated signal with ao = 0.5, and a positively correlated 
signal with ao = 0.8. We find that in the intermediate 
regime F{n) (obtained from all methods) exhibits an ap- 
proximate power-law behavior characterized by a single 
scaling exponent a. At large scales n for DFA-0, DMA, 
and WDMA-1, we observe a crossover in F{n) leading to 
a flat regime. With increasing ao this crossover becomes 
more pronounced and moves to the intermediate scal- 
ing range. In contrast, such a crossover at large scales 
is not observed for DFA-1, indicating that the DFA-1 
method can better quantify the correlation properties at 
large scales. At small scales n the scaling curves F{n) ob- 
tained from all methods exhibit a crossover which is more 
pronounced for anti-correlated signals (ao — 0.2) and be- 
comes less pronounced for uncorrelated (ao = 0.5) and 
positively correlated signals (ao = 0.8). 

We next systematically examine the performance of 
the DFA-0, DFA-1, DMA and WDMA-1 methods by 
varying ao over a very broad range of values (0.1 < ao < 
3.5) [Fig. 12]. For all four methods, we compare ao with 
the exponent a obtained from fitting the rms fluctuation 
function F{n) in the scaling range 10^ < n < 10^, i.e., 
the range where all methods perform well according to 
our observations in Fig. ^ If the methods work prop- 
erly, for each value of the "input" exponent ao we expect 
the estimated "output" exponent to be a — ao. We 
find that the scaling exponent a, obtained from different 
methods, saturates as the "input" correlation exponent 
ao increases, indicating the limitation of each method. 
The saturation of scaling exponent at a = 1 indicates 
that DMA and WDMA-1 do not accurately quantify the 
correlation properties of signals with ao > 1. 

In contrast, the DFA-^ method can quantify accurately 
the scaling behavior of strongly correlated signals if the 
appropriate order £ of the polynomial fit is used in the 
detrending procedure. Specifically, we find that the val- 
ues of the scaling exponent a obtained from the DFA-£ 
are limited to a < £ + 1. Thus the DFA-£ can quan- 
tify the correlation properties of signals characterized by 
exponent ao < £ + I. For signals with ao > ^ -I- 1 we 
find that the output exponent a from the DFA-^ method 
remains constant at a = £ + 1. These findings suggest 
that in order to obtain a reliable estimate of the correla- 
tions in a signal one has to apply the DFA-£ for several 
increasing orders £ until the obtained scaling exponent a 
stops changing with increasing £. 

Since the accuracy of the scaling exponent obtained 
from the different methods depends on the range of scales 
n over which we fit the rms fiuctuation function F(n) (as 
seen in Fig.^, and since different methods exhibit differ- 
ent limitations for the range of scaling exponent values 
(as demonstrated in Fig. [^J , we next investigate the local 
scaling behavior of the F(n) curves to quantify the per- 
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FIG. 3: A comparison of the local scaling exponent aioc as a function of the scale n for the DMA, WDMA-1, DFA-0, and 
DFA-1 methods. We consider signals of length = 2'^^ and varying values of the correlation exponent ao. The local scaling 
exponent aioc quantifies the stability of the scaling curves F{n) (see Fig.0 and is expected to exhibit small fluctuations around 
a constant value ao if F{n) is well fitted by a power-law function, qq is denoted by horizontal dotted lines. Symbols denote 
the estimated values of aioc and represent average results from 50 realizations of artificial signals for each value of the "input" 
scaling exponent ao. Deviations from the horizontal lines at small or at large scales indicate limitations of the methods to 
accurately quantify the build-in correlations in different scaling ranges. 



formancc of the different methods in greater details. To 
ensure a good estimation of the local scaling behavior, we 
calculate F{n) at scales n = 4x2*/^"', i = 0,1, 2, which 
in log scale provides 64 equidistant points for F{n) per 
bin of size log 2. To estimate the local scaling exponent 
aioc, we locally fit F{n) in a window of size w = 3 log 2, 
e.g, aioc is the slope of F{n) in a window containing 3 x 64 
points. To quantify the detailed features of the scaling 
curve F(n) at different scales n, we slide the window w 
in small steps of size ^ — \ log 2 starting at n = 4, thus 
obtaining approximately 70 equidistant aioc in log scale 
per each scaling curve. We consider the average value 
of aioc obtained from 50 different realizations of signals 
with the same correlation exponent ao. 

In Fig. 121 we compare the behavior of aioc as a function 
of the scale n to more accurately determine the best fit- 
ting range in the scaling curves F{n) obtained from the 
DMA, WDMA-1, DFA-0, and DFA-1. A rnis fluctua- 



tion function exhibiting a perfect scaling behavior would 
be characterized by aioc — ao for all scales n and for 
all values of ao denoted by horizontal lines in Fig. O A 
deviation of the aioc curves from these horizontal lines in- 
dicates an inaccuracy in quantifying the correlation prop- 
erties of a signal and the limitation of the methods. Our 
results show that the performance of different methods 
depends on the "input" ao and scale n. At small scales 
and for ao < 0.8 we observe that aioc for all methods 
deviates up from the horizontal lines suggesting an over- 
estimation of the real correlation exponent ao. This ef- 
fect is less pronounced for uncorrelated and positively 
correlated signals. At intermediate scales aioc exhibits 
a horizontal plateau indicating that all methods closely 
reproduce the input exponent for ao < 0.8. This in- 
termediate scaling regime changes for different types of 
correlations and for different methods. At large scales 
of n > 10^, the DMA, WDMA-1, and DFA-0 methods 
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FIG. 4: Values of the local scaling exponent aioc as a function 
of the scale n obtained from 20 different realizations of artifi- 
cial anti-correlated signals with an identical scaling exponent 
ao = 0.2. 



FIG. 6: Values of the local scaling exponent ai^c as a func- 
tion of the scale n obtained from 20 different realizations of 
artificial positively correlated signals with an identical scaling 
exponent ao = 0.8. 
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FIG. 5: Values of the local scaling exponent aioc as a function 
of the scale n obtained from 20 different realizations of arti- 
ficial uncorrelated signals with an identical scaling exponent 
ao = 0.5. 



strongly underestimate the actual correlations in the sig- 
nal, with aioc curves sharply dropping for all values of 
ao [Fig. EIa),(b),(c)]. In contrast, the DFA-1 method 
accurately reproduces ao at large scales with aioc fol- 
lowing the horizontal lines up to approximately TV/ 10 
[Fig- Eld)]. In addition, the DFA-1 method accurately 
reproduces the correlation exponent at small and inter- 
mediate scales even when ag > 1 [Fig. Old)], while the 
DMA, WDMA-1 and DFA-0 are limited to ao < 0.8. 

For a certain "input" correlation exponent ao, we can 
estimate the good fitting regime of F{n) to be the length 
of the plateau in Fig. 13 For example, for ao — 0.2 the 
calculated scaling exponent aioc obtained from the DMA 
method is approximately equal to the expected value 
ao = 0.2 within a range of two decades (10^ < n < 10'*). 
Similarly, the good fitting range of F{n) obtained from 
the DFA-0 for ao = 0.2 is about three decades (10^ < 
n < 10^). However, the calculated local scaling expo- 
nent aioc can fluctuate for different realizations of cor- 
related signals. Although the mean value obtained from 
many independent realizations is close to the expected 
value, the fluctuation of the estimated scaling exponent 
can be very large. Thus, it is possible for aioc to deviate 
from ao and the scaling range estimated from Fig.|2|may 
not be a good fitting range. Therefore, it is necessary 
to study the dispersion of the local scaling exponent to 
determine the reliability of the "good" fitting range esti- 
mated from Fig.O In Fig.^|Slin|we show the results for 
aioc from 20 different realizations of the correlated signal 



7 



10' 



10' 



10' 



10 



(a) DMA 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



10' 



10 



10' 



10' 



(c)WDMA-1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



(b) DFA-0 




I 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 





1 


1 

0.9 


10= 






0.8 








0.7 


10^ 






0.6 








0.5 








0.4 


10= 






0.3 








0.2 


10^ 






0.1 














(d) DFA-1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



I 



FIG. 7: Probability density of the estimated values of ao — ^ < aioc < ao + S, where S = 0.02 for a varying scale range n and for 
different values of the "input" correlation exponent oq. Separate panels show the performance of the DMA, WDMA-1, DFA-0 
and DFA-1 methods, respectively, based on 50 realizations of correlated signals for each value of ao- The probability density 
values p are presented in color, with the darker color corresponding to higher values as indicated in the vertical columns next 
to each panel. A perfect scaling behavior would correspond to dark-colored columns spanning all scales n for each value of uq. 



with ao = 0.2, ao = 0.5, and ao = 0.8 respectively. For 
all methods, we observe that there is a large dispersion of 
aioc, indicating strong fluctuations in the scaling function 
F{n) at large scales n (n - 10^ for DMA and WDMA-1 
and n - 10'' for DFA-0 and DFA-1) [Fig. H EJ 0. This 
suggests that the good fitting range obtained only from 
the mean value of aioc, as shown in Fig. |2| may be over- 
estimated. 

To better quantify the best fitting range for different 
methods and for different types of correlations we de- 
velop a three-dimensional representation [Fig.[7]- Based 
on 50 realizations of correlated signals with different val- 
ues of 0.1 < ao < 1.1, for each scale n we define the 
probability p (normalized frequency) to obtain values for 
ao — 6 < aioc < ao + (5, where S — 0.02 (arguments sup- 
porting this choice of 5 are presented in section 111(B)). 
Again, as in Fig.|31 for each realization of correlated sig- 
nals with a given ao, we calculate aioc by fitting the rms 
fluctuation function F{n) in a window of size w = 3 log 2 
sliding in steps of A = ■|log2. Vertical color bars in 
Fig. represent the value of the probability p — darker 



colors corresponding to higher probability to obtain ac- 
curate values for aioc- Thus dark-colored columns in the 
panels of Fig. [7| represent the range of scales n where the 
methods perform best. 

For the DMA and WDMA-1 methods, we find that 
with high probability {p > 0.7), accurate scaling results 
can be obtained in the scaling range of two decades for 
0.4 < ao < 0.6. However, WDMA-1 performs better 
at small scales compared to DMA. For explanation why 
the WDMA-1 performs better at small scales compared 
to DMA, see Appendix II. In contrast, DFA-0 exhibits 
an increased fitting range of about three decades for 
0.4 < ao < 0.8, while for the DFA-1 we find the best fit- 
ting range to be around three decades for ao > 0.5. For 
strongly anti-correlated signals (ao < 0.2), all methods 
do not provide an accurate estimate of the scaling expo- 
nents ao. However, by integrating anti-correlated signals 
with ao < 0.3 and applying the DFA-1 method, we can 
reliably quantify the scaling exponent, since DFA-1 has 
the advantage to quantify signals with ao > 1 [Fig.^d)]. 
This can not be obtained by the other three methods 
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[Fig.|7|;a),(b),(c)]. 



B. Centered moving average method and DFA 

To test the accuracy of the CDMA method we per- 
form the same procedure as shown in Fig|2| We calculate 
the local scaling exponent ajoc for signals with different 
"input" correlation exponent ao and for a broad range 
of scales n [Fig|H]. We find that for 0.3 < ao < 0.8 
the CDMA method performs better than the DMA for 
all scales n, and the average value of aioc follows very 
closely the expected values of a indicated by horizontal 
lines in Fig|Sl For anti-correlated signals with ao < 0.3, 
both DMA and CDMA overestimate the value of ao 
at small scales n < 10^. For strongly correlated sig- 
nals with ao > 0.8, CDMA underestimates ao at small 
scales n < 10^, in contrast to DMA which overestimates 
ao. For correlated signals with ao > 1.1 (not shown in 
FiglSJc)) the deviation of ajoc from the expected value ao 
for the CDMA method becomes even more pronounced 
and spreads to large scales. At intermediate and large 
scales CDMA performs much better — aioc closely fol- 
lows the horizontal lines [Fig|Sl(a),(c)]. These differences 
in the performance of the DMA and CDMA methods are 
also clearly seen in the probability density plots shown 
in Figli 

Next, we compare the stability of the DMA, CDMA, 
DFA-0, and DFA-1 methods in reproducing the same "in- 
put" value of ao for different realizations of correlated 
signals. We generate 50 realizations of signals for each 
ao, and we obtain the average value and the standard de- 
viation of aioc for a range of scales n. The values of the 
standard deviation are represented by error bars in Fig|Sl 
for each value of aioc at all scales n. We find that with 
increasing scales n, the standard deviation gradually in- 
creases, and that for DMA the standard deviation is less 
than 0.02 while for DFA the standard deviation is less 
than 0.01 in the range of scales n up to iV/100 {N is the 
signal length). For all methods at scales n > iV/100, the 
standard deviation increases more rapidly, and thus the 
stability of the methods in reproducing the same value 
of the exponent for different realizations decreases. 

In Fig^l we present the dependence of aioc on the 
scale n for the weighted centered detrended moving av- 
erage method. Compared to the CDMA method, the 
WCDMA method weakens the overestimation of a/oc at 
small scale for anti-correlated signals and provides accu- 
rate results of a/oc at small scales for positively corre- 
lated signals with 0.5 < ao < 1. Compared to the DFA 
method, the WCDMA performs better at small scales 
for 0.5 < ao < 1.0. However, at larger scales n > 10^, 
the standard deviation of DFA-1 is smaller than that of 
WCDMA (Figs|HId),[iniand[ni), indicating more rehable 
results for the local scaling exponent aioc obtained from 
DFA-1. 

Finally, we test how the choice of the parameter S will 
affect the probability density plots shown in Fig[7| and 



FigEl To access the precision of the methods one has to 
increase the confidence level by decreasing S. In Fig[7|and 
Figl^lwe have chosen 6 = 0.02 to correspond to the value 
of the standard deviation for aioc at scales n < 10* as 
estimated by the DMA method [Fig|S]. We demonstrate 
that the distribution plot for DMA with S = 0.02 (shown 
in Fig[7|l changes dramatically when we chose S = 0.01 
(as shown in Fig ll2r b)). This result confirms the obser- 
vation from FiglS^a) and (d) that the DFA-1 method is 
more stable (smaller standard deviation) and more accu- 
rate (average of a/oc closer to the theoretically expected 
value ao) than the DMA method. 



IV. DISCUSSION 

We have systematically studied the performance of the 
different variants of DMA method when applied to sig- 
nals with long-range power-law correlations, and we have 
compared them to the DFA method. Specially, we have 
considered two categories of detrended moving average 
methods — the simple moving average and the weighted 
moving average — in order to investigate the effect of 
the relative contribution of data points included in the 
moving average window when estimating correlations in 
signals. To investigate the role of "past" and "future" 
data points in the dynamic averaging process for signals 
with different correlations, we have also considered the 
cases of backward and centered moving average within 
each of the above two categories. Finally, we have in- 
troduced a three-dimensional representation to compare 
the performance of different variants of the DMA method 
and the DFA methods over different scaling ranges based 
on an ensemble of multiple signal realizations. 

We find that the simple backward moving average 
DMA method and the weighted backward moving av- 
erage method WDMA-^ have limitations when applied 
to signals with very strong correlations characterized by 
scaling exponent ao > 0.8. A similar limitation is also 
found for the £ = order of the DFA method. How- 
ever, for higher order £, the DFA-£ method can accu- 
rately quantify correlations with ao < £ + 1. We also 
find that at large scales the DMA, WDMA-£, and DFA-0 
methods underestimate the correlations in signals with 
0.5 < ao < 1.0, while the DFA-^ method can more ac- 
curately quantify the scaling behavior of such signals. 
Further, we find that the scaling curves obtained from 
the DFA-1 method are stable over a much broader range 
of scales compared to the DMA, WDMA-1, and DFA- 
methods, indicating a better fitting range to quantify 
the correlation exponent ao. In contrast, we find that 
WDMA-£ with a higher order £, more accurately repro- 
duce the correlation properties of anti-correlated signals 
(ao < 0.5) at small scales. Accurate results for anti- 
correlated signals can also be obtained from the DFA-1 
method after first integrating the signal and thus reduc- 
ing the value of the estimated correlation exponent by 1. 

In contrast to the simple backward moving average 
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FIG. 8: A comparison of the local scaling exponent aioc as a function of the scale n for the DMA, CDMA, DFA-0, and DFA-1 
methods. We consider signals of length A*' — 2^" and varying values of the correlation exponent ao. The local scaling exponent 
oioc quantifies the stability of the scaling curves F{n) and is expected to exhibit small fluctuations around a constant value ao 
if F{n) is well fitted by a power-law function, ao is denoted by horizontal dotted lines. Symbols denote the estimated values 
of Qioc and represent average results from 50 realizations of artificial signals for each value of the "input" scaling exponent qq. 
Deviations from the horizontal lines at small or at large scales indicate limitations of the methods to accurately quantifying 
the build-in correlations in different scaling ranges. Error bars represent the standard deviation for each average value of aioc 
at different scales n, and determine the accuracy of each method. 



(DMA) and DFA-0 methods, the centered moving av- 
erage CDMA provides a more accurate estimate of the 
correlations in signals with 0.3 < ao < 0.7 at small 
scales n < 10^, and in signals with ao > 0.7 at inter- 
mediate scales 10^ < n < 10^. However, the CDMA 
method strongly underestimates correlations in signals 
with ao > 0.7 at small scales (n < 10^), while the DFA-1 
method reproduces quite accurately the correlations of 
signals with ao > 0.7 at both small and intermediate 
scales. We also find that by introducing weighted cen- 
tered moving average WCDMA, one can overcome the 
limitation of the CDMA method in estimating correla- 
tions in signals with ao > 0.5 at small scales {n < 10^). 
On the other hand, the WCDMA method is characterized 
by larger error bars for a/oc at intermediate scales com- 
pared to the CDMA method. Further, we find that the 
performance of the WCDMA is comparable to the DFA- 
1 method for signals with 0.5 < ao < 1. At small scales 



the WCDMA performs better than the DFA-1 method, 
while at the intermediate scales 10^ < n < 10'', DFA-1 
provides more reliable local scaling exponent with smaller 
standard deviation based on 50 independent realizations 
for each ao. For very strongly correlated signals with 
ao > 1, we find that the DFA-1 method performs much 
better at all scales compared to WCDMA and all other 
variants of the DMA method. 
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FIG. 9: Probability density of the estimated values of ao — S < 
Qfioc < ao + S, where 5 = 0.02 for a varying scale range n and 
for different values of the "input" correlation exponent ao- 
The two panels show the performance of the DMA and CDMA 
methods, respectively, based on 50 realizations of correlated 
signals for each value of ao- The probability density values 
•p are presented in color, with the darker color corresponding 
to higher values, as indicated in the vertical columns next to 
each panel. A perfect scaling behavior would correspond to 
dark-colored columns spanning all scales n for each value of 



V. APPENDIX I 

Higer order weighted moving average 

To account for different types of correlations in sig- 
nals, we consider the ^-order weighted moving average 
(WDMA-^), defined as 

y„W = ^^^y(z-fc) + Ay„(z-^), (13) 

fe=0 

where X — nj (n-\-l\ i is the order of the moving average, 
y(i) is defined in Eq.(Q. The relative importance of the 
two terms entering the function in Ea. (|13|l . can be further 
understood by analyzing the properties of the transfer 
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n 

FIG. 10: Local scaling exponent oioc as a function of the scale 
n for the WCDMA method. We consider signals of length 
N = 2^° and varying values of the correlation exponent ao- 
The expected value of the exponent ao is denoted by hori- 
zontal dotted lines. Symbols denote the estimated values of 
aioc and represent average results from 50 realizations of ar- 
tificial signals for each value of the "input" scaling exponent 
ao. Deviations from the horizontal lines at small or at large 
scales indicate limitations of the methods to accurately quan- 
tify the build-in correlations in different scaling ranges. Error 
bars represent the standard deviation for each average value 
of ctioc at different scales n, and determine the accuracy of 
the method. 



function H{f) in the frequency domain (see Appendix 
II). 

Compared to the traditional exponentially weighted 
moving average (of order i = 1 ) where the terms in 
Ea. (|13|) decrease exponentially, the higher order i > 1 
allows for a slower, step-size decrease of the terms in 
Ea. (|13() with a "step" of size £. The fluctuation func- 
tion F{n) is obtained following Ea. ((Tn|l and Ea. ((TT|l . The 
WDMA-^ allows for a more gradual decrease in the dis- 
tribution of weights in the moving average box, and thus 
may be more appropriate when estimating the scaling 
behavior of anti-correlated and uncorrelated signals. 

We apply the WDMA-f method for increasing values 
of i to correlated signals with varied values of the scaling 
exponent ao- To study the performance of the WDMA- 
£ methods, we estimate the scaling behavior of the rms 
fluctuation function F{n) at different scales n by calcu- 
lating the local scaling exponent aioc in the same way 
as discussed in Fig. O We find that at large scales for 
£ = 2, 5, the aioc curves deviate significantly from the 
expected values ao — presented with horizontal dashed 
lines in Fig.[T21 This indicates that the WDMA-^ method 
significantly underestimates the strength of the correla- 
tions in our artificially generated signals. Further, as for 
£ = 1, we find that for higher order £ > 1 the WDMA- 
£ methods exhibit an inherent limitation to accurately 
quantify the scaling behavior of positively correlated sig- 
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FIG. 11: A comparison of the local scaling exponent aioc 
as a function of the scale n obtained from (a) the WCDMA 
method and (b) the DFA-1 method. Symbols denote the esti- 
mated values of aioc calculated as in Figlinifor different " m- 
put" scaling exponents ao > 1. Error bars representing the 
standard deviation around the average aioc are smaller for 
the DFA-1 method at all scales n, indicating that the DFA-1 
method provides more reliable results. 

nals with ao > 0.7. This behavior is also clear from 
our three-dimensional presentation in Fig. 1141 For anti- 
correlated signals, however, the WDMA-£ performs bet- 
ter at small and intermediate scales for increasing order 
i as ao decreases [Fig. (see Appendix II). These ob- 
servations are also confirmed from the three-dimensional 
probability histograms in Fig. 1141 where it is clear that 
the scaling range for the best fit shrinks for positively 
correlated signals (ao > 0.5) for increasing order £, while 
for anti-correlated signals (ao < 0.5), there is a broader 
range of scales over which a best fit (with a probability 
of p > 0.7) is observed. 

VI. APPENDIX II 

Moving average methods in frequency domain 

In this appendix, the performance of the DMA algo- 
rithm is discussed in the frequency domain. The interest 
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FIG. 12: Probability density of estimated values of ao — 5 < 
Qfioc < ao + S, where 5 = 0.01 for varying scale range n and 
for different values of the "input" correlation exponent qq. 
The two panels show the performance of the DMA and DFA- 
1 methods, respectively, based on 50 realizations of correlated 
signals for each value of ao. The probability density values 
p are presented in color, with darker color corresponding to 
higher values as indicated in the vertical columns next to each 
panel. A perfect scaling behavior would correspond to dark- 
colored columns spanning all scales n for each value of qq. 



of the frequency domain derives from the simplification 
designed to describe the effect of the detrending function 
yn{i) in terms of the product of the square modulus of the 
transfer function ff„(/) and of 5'(/), the power spectral 
density of the noisy signal y{i). 

The simple moving average j/n(?) of window size n is 
defined as 

^ n — r 

n ^-^ 

k=0 

corresponding to the discrete form of the causal convolu- 
tion integral, where the convolution kernel introduces the 
memory effect. Ea. p4f) is a sum with a constant mem- 
ory kernel h{t), i.e., a step function with an amplitude 
1/n [Fig ll5f a'l]. The function h{t) uniformly weights the 
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FIG. 13; A comparison of the local scaling exponent aioc as a function of the scale n for the WDMA-^ method with different 
order £ — 2, 5 of the weighted moving average. We consider signals of length A*" — 2^" and varying values of the correlation 
exponent qq. The local scaling exponent aioc quantifies the stability of the scaling curves F{n) (see Fig.0, and is expected to 
exhibit small fluctuations around a constant value qq if F{n) is well fitted by a power-law function, qo is denoted by horizontal 
dotted lines. Symbols denote the estimated values of Qioc and represent average results from 50 realizations of artificial signals 
for each value of the "input" scaling exponent uq. For small values of i at small and intermediate scales n, WDMA-^ accurately 
reproduces the scaling behavior of signals with 0.4 < ao < 0.8, while for large £ , the scaling behavior of anti-correlated signals 
with ao < 0.4 are better reproduced at small scales. 



contribution of all the past events in the window [0, n), 
thus it works better for random paths with a correla- 
tion exponent centered around 0.5. For higher degrees of 
correlation/anti-correlation, it should be taken into ac- 
count (as already explained in the section describing the 
DMA function) that each data point is more correlated 
to the most recent points than to the points further away. 

In the frequency domain, y„{i) is characterized by the 
transfer function Hn{f) (the Dirichlet kernel), which is 

F„(/) = ^^-e— /. (15) 
nnf 

Hn{f) takes the values i?„(0) = 1 and i?„(fc/o) = for 
k = l,2, ...n. 

The transfer function H{f) of any filter should ideally 
be a window of constant amplitude, going to zero very 
quickly above the cut-off frequency l/n. By observing 
the curves of Fig llSf b') and Fig ll(jl it is clear that the 



filtering performance of Hn (/) is affected by the presence 
of the side lobes at frequency higher than l/n. 

As observed in Fig^| |i/4(/)|^ presents a side lobe 
allowing the components of the signal y{i), with a fre- 
quency between l/n and 2/n (i.e., time scales between 
n/2 and n) to pass through the filter, thus giving a spuri- 
ous contribution to y„(i). These components contribute 
to the rms F{n) (defined in Ea. (|ll|l ) less than what 
should correspond to n on the basis of the scaling law 
F{n) ~ , with the consequence of increasing the slope 
F{n) at small scales. 

We next discuss the reasons why the weighted mov- 
ing average might reduce this effect. The exponentially 
weighted moving average (WDMA-^) weights recent data 
more than older data. It is defined as 

y„,,(i) = ■^i^^y(i-fc) + Ay„,,(i-^) . (16) 

fe=0 
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(d) WDMA-5 
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FIG. 14: Probability density of estimated values of oq — (5 < aioc < ao + 5, where 5 — 0.02 for the varying scale range 
n and for different values of the "input" correlation exponent ao- Separate panels show the performance of the WDMA-2, 
WDMA-3, WDMA-4, and WDMA-5 methods, respectively, based on 50 realizations of correlated signals for each value of ao. 
The probability density values p are presented in color, with the darker color corresponding to higher values, as indicated in 
the vertical columns next to each panel. A perfect scaling behavior would correspond to dark-colored columns spanning all 
scales n for each value of ao. 



The coefficients are commonly indicated as weights of the 
filter and are given by 



A = 



(17) 



Taking the Fourier transform on Eq. 116(1 . we obtain 

YnAf) = (1 - • Yif) + \%j{f) ■ e-'^^'f, (18) 

where Y{f), Yn^f) are the Fourier transforms of y{i) 
and yn/ii) respectively. Further, we have 



YnAf) 



1-A 



1 - Ae-»2^^/ 

Thus the transfer function is 

1 - A 



H,{f)-Y{f). (19) 



HnAf) 



1 - Ae-»27rf 



Hiif)- 



(20) 



From Ea. (|2U|l . one can find that the cutoff frequency for 
\HnAf)\^ is min[l/{2ny/n{n + £)),!/£]. In Fig.[li the 
transfer function of the weighted moving averages with 
n = 4 and £ = 2 and with n = 4 and i = A respectively 
are shown. It can be observed that the effect of the side 
lobe to the performance of |i?4,2(/)|^ and |_ff4.4(/)|2 has 
become negligible compared to that of |iJ4(/)p , with the 
consequence of reducing the high frequency components 
in the detrended signal and thus reducing the deviation 
of the Qfioc, as discussed in the paper. 
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FIG. 16: Plot of the function \H{f)\^ for the simple mov- 
ing average with n = 4, \H4{f)f; for the weighted moving 
average, with n — 4 and £ = 2, |^^4,2(/)P respectively ; for 
the weighted moving average with n = 4 and £ — 4 1 1/4,4 
respectively. 
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